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Abstract. Recent data, including the three-year WMAP data, the full 2dF galaxy 
power spectrum and the first-year data of the Supernova Legacy Survey, are used to 
constrain model parameters in quintessence cosmologies. In particular, we discuss the 
inverse power-law (RP) and SUGRA potentials and compare parameter constraints 
with those for ACDM. Both potentials fit current observations with a goodness of fit 
comparable or better than ACDM. The constraints on the energy scale Ade appearing 
in both potential expressions are however different. For RP, only energy scales around 
the cosmological constant limit are allowed, making the allowed models quite similar 
to ACDM. For SUGRA, Ade values approximately up to Electroweak energy scale are 
still allowed, while other parameter intervals are slightly but significantly displaced. 
In particular a value of the primeval spectral index n s = 1 is still allowed at the 95 % 
c.l., and this can have an impact on constraints on possible inflationary potentials. 
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1. Introduction 

In the last two decades, observations of high redshift supernovae pQ showed that the 
expansion of the Universe is undergoing a phase of positive acceleration. Together with 
measurements of Large Scale Structure |2j and of anisotropies of the Cosmic Microwave 
Background (CMB) temperature and polarization these data can be explained if 
about 70 % of the energy density of our universe is made up by a component with 
negative pressure, dubbed Dark Energy (DE). 

The basic candidate for DE is the cosmological constant, A, with an equation of 
state w = p/ p = — 1. As showed by the three- year Wilkinson Microwave Anisotropy 
Probe (WMAP) observations jlj, a power-law ACDM model with adiabatic initial 
conditions is a good fit to most cosmological data [3] . Within the framework of a ACDM 
cosmology, the combination of the three-year WMAP (WMAP3) observations with LSS 
and SNIa data, allowed to measure several fundamental cosmological parameters with 
percent accuracy. However, a cosmological constant has severe problems on theoretical 
grounds, in particular fine-tuning and cosmic coincidence. 

Several alternatives to a cosmological constants have been proposed, including a 
perfect fluid with constant self-interacting scalar fields, with a time depending 

w > — 1 (quintessence [3 El) or w < — 1 (phantom [Hj), or with non-canonical kinetic 
term (k-essence [Hj; quintom [TUj): modifications of General Relativity [11J. 

Quintessence models are particularly interesting due to the existence of tracker 
potentials. Under the proper assumptions for the potential [12], the energy density of 
quintessence follows closely the evolution of the dominant component, though decreasing 
more slowly. In addition, the scalar field evolution does not depend on the initial 
conditions, lessening fine-tuning and related problems. On the contrary, dynamics 
depend on the current energy density of DE and on the shape of the scalar field potential 
V((f>) and its parameters, which should be added to the usual cosmological ones during 
data analysis. 

Several authors used a combination of WMAP3 data and other cosmological data, 
to study DE properties. The WMAP team placed constraints on the value of w, 
assumed constant, while other works studied the evolution of pde [HI]- or assumed 
phenomenological parameterization for w(a) [EJ- However, the resulting constraints 
may depend on the choice of parameterization (e.g. an d, i n general, do not translate 
directly into constraints on the shape of the potential. 

In this work we use recent cosmological data to directly constrain the shape of 
two classes of quintessence potentials: the Ratra-Peebles Potential (RP ^Hj) and the 
SUGRA potential [Zj . The plane of the paper is as follow. In Sec. El we outline the 
models studied together with the data and the method used. In Sec. El we describe our 
results and discuss them in Sec. HI 
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2. Models and Data 

In the limit of spatial homogeneity, the equation of motion for a quintessence scalar field 
<p is given by 

+ m'<\> + a 2 V\<f>) = (1) 

where dots denote derivatives with respect to conformal time r, the Hubble parameter 
H = a I a, and we used natural units, c = h = 1. The energy and pressure associated 
with the field are given by 





+ V{(t>) ; P4t = IZ \ - V {(j>) . (2) 



Together with Friedman equations, (£[]) and (0) allow to solve for the background 
evolution in a quintessence cosmology, once the the potential V((p), and the density 
parameters of the different components fib, fic, etc., are assigned. At variance with a 
cosmological constant, quintessence is not a smooth component and perturbations in 
the scalar field need to be accounted for [17J. In the Synchronous gauge they obey the 
equation 

+ 2H<p + a 2 V"{<f>) + ^ + V i (p = O , (3) 

where <p = 5(f) and f) is the trace part of the metric perturbation. Perturbation evolution 
is also fixed by the potential. 

Quintessence cosmologies, then, clearly depend on shape of V((p), which in principle 
is totally arbitrary. Considering only tracker potentials allow to ease the fine tuning 
problem of ACDM cosmologies. As said in Sec. ^ in this work we focus on two classes 
of tracker potentials, RP and SUGRA: 

V{<t>) = RP (4) 



nt>) = ^ exp (^-j , SUGRA (5) 

here m p is the Planck mass. The former was one of the first tracker potentials studied, 
and is characterized by a slowly varying equation of state. The latter originates within 
the context of Supergravity theories and the corresponding equation of state undergoes 
a rapid transition when the field becomes dominant. Both potentials depend on two 
parameters, the slope a and the energy scale A^e- For a given choice of parameters, the 
scalar field associated with the two potentials evolve in a similar way until the epoch 
of DE domination; at that point in a SUGRA model w(a) rapidly falls toward values 
— 1 < w(l) < —0.9, while in a RP model w(a) decreases more slowly, and stays at 
considerably higher values. 

In a cosmological context, a and Ade are not independent; fixing the current energy 
density of DE, pde,o, leaves only one free parameter. Notice that for a ~ 0, we recover 
the behaviour of the cosmological constant. Here, we take A^e as a free parameter; 
however, different conventions are possible fHl OS] ■ 
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In this work, we test RP and SUGRA models using recent data, and try to place 
constraints on the value of A^e- We assume a power-law spectrum of primordial 
density fluctuations with adiabatic initial conditions, and restrict our analysis to flat 
cosmologies with no massive neutrinos, so that ^de — 1 — £\ — fi c - The models 
studied are thus specified by seven parameters: the physical baryon, iO\> = ^h 2 , and 
CDM, u c , densities; the angular size of the sound horizon at recombination 6 S ; the 
slope, n s , and amplitude, A s , of the power-spectrum of density fluctuations; the optical 
depth to reionization, r; and A = Log 10 (ADE/GeV). Notice that h is the reduced 
Hubble parameter, H = 100/iKm/s/Mpc. This choice of parameters is particularly 
efficient in minimizing the impact of the degeneracies of CMB angular power spectra 
on data analysis [20] • To explore the parameter space, we used a modified version of 
the COSMOMC package. We assumed a flat priors on all parameters; in particular 
—11.75 < A < 19. The upper limit corresponds to the Planck mass, while the lower 
one is set by (A DE ) 4 ~ lCT 47 GeV 4 , i.e. the cosmological constant energy density in a 
concordance ACDM model. We also tested that a different choice of priors, in particular 
using H instead of 9 S , does not significantly alter our results. 

Models were tested against WMAP3 data, the galaxy power spectrum P(k) of the 
full 2dF Galaxy Redshift Survey i2l. and the distance measurements to high-redshift 
supernovae by the Supernova Legacy Survey (SNSL 22]). Assuming the three sets of 
data to be independent, the total likelihood is the product of the individual likelihoods 

£tot = £ C mb x £ LSS x £ S n • (6) 

We also define an effective \ 2 = — 21n(£). 

In the analysis of the CMB anisotropy spectrum, the WMAP team included the 
contribution of secondary anisotropies due to Sunayev-Zel'dovich (SZ) effect 0. The 
SZ contribution was parametrized in term of the amplitude of the signal, Asz, relative 
to a reference ACDM model and final constraints were obtained after marginalization 
over Asz- The effects of SZ marginalization on ACDM estimates are very small. In 
addition, the SZ contribution depends on the mass function and its evolution; these, 
in turn, depend on the quintessence model considered and require extensive numerical 
simulations [23] . Thus, we set Asz = here. We also did not include small scale data in 
our analysis of CMB anisotropies. High-£'s mostly probe the recombination history and 
the running of the spectral index, but are not significantly affected by the quintessence 
models we considered. 

Analysis of LSS data requires the introduction of an additional parameter, b, 
accounting for the bias in the clustering of galaxies with respect to the clustering 
of the matter field. Here, we assume that b is independent of scale in the range of 
0.02 /iMpc -1 < k < 0.15 /iMpc -1 considered [21], and the code performs an analytic 
marginalization over it. In these range of wavenumbers non-linear corrections are 
small [21], and we do not include them in the analysis. For both SUGRA and RP 
potentials we run 10 chains of ~ 50000 points each. We also run a set of chains for a 
ACDM model. Comparing results for quintessence with those for ACDM allows then to 
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Table 1. SUGRA parameters. For all parameters except Ade we show the mean value 
(x) = J dx'£(x')x' and corresponding 68% confidence level interval, for the different 
combinations of datasets used. For Ade we show the 95% upper confidence limit. 



Parameter 


WMAP3 


WMAP3 + 2dF 


WMAP3 + SNLS 


WMAP3 + 2dF 
+SNLS 


100Q b h 2 


o ocro+0.077 
Z - ZJZ -0.078 


2.2461K 


9 OQ/1 +0.074 
Z.ZO^_oo79 


o 9/1 9+0.073 
z - z ^ z -0.077 


Q c h 2 


n 1 n^o+o- 0082 

U.±UOZ_ .oo82 


0.09801°;^ 


1001 +aoon 

U.1UU1_ 0075 


0.09911°;^ 


Ho 




67.6±l;f 


70.6i;i 


70.1111 


T 

n s 


093 +a014 

967 +a020 
u.yu i _ .020 


098 +a015 

U.UJO_ 015 

96 c i +a019 
u.yuo_ 019 


097+ ' 015 
u.uy i _ .oi5 

u.yoi_ .oi8 


098+ - 015 
u.uyo_ 015 

963+ 017 


Log 10 (A DE ) 


< 14.9 


< 9.3 


< 0.1 


< 2.1 




n 09 e: +0.070 


0.26618S 


n 247 +0.024 
u - z ^' -0.024 


0-248±g : gJ 


^8 


n fiS4 +0 085 
u.uo<±_ a086 


o.648±g:8S 


696+ - 065 

U.UJU_o.063 


0.685±g : ggg 



factor out possible biases introduced by our choice of datasets, priors and methods. 

3. Results 

Table Q and table 121 summarizes our main results for SUGRA and RP cosmologies, 
respectively. For each parameter x, except A, we list the mean value 

(x) = J dx'C{x')x' (7) 

and the limits of the 68% confidence interval, defined as appropriate quantiles of the 
marginalized distribution, for different combinations of data sets. For A we list the 
upper 95% confidence limit. Table El instead, shows results for a basic 6-parameters 
ACDM cosmology, in order to outline the effects of quintessence. 

In figure Qand figure H3 we plot the marginalized ID likelihoods, for some basic and 
derived parameters. Different colours and lines types refer to different combinations of 
data (see caption for details). In figure El and figure HI we show the joint 2D confidence 
limits between A DE and other relevant parameters. 

4. Discussion 

Using the first-year WMAP (WMAP1) data, several authors ^H] placed constraints on 
the slope of RP potential, finding a < 2. This constraint does not directly translate into 
a constraint on the energy scale A DE , as the a-A DE relation depends on Pde,o- However, 
it was clear that allowed values fall well below the 1 — 10 6 GeV range, which includes 
the Electroweak (EW) and (possibly) the Supersymmetry (SUSY) breaking scales. Let 
us outline soon that WMAP3 data do not improve this situation. While a SUGRA 
potential still provides a fair fit of data for a physically relevant range of Ad E values, 
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Figure 1. Marginalized ID constraints on basic and derived parameters of 
SUGRA cosmology. Dotted lines show constrains from WMAP3 only, dashed lines 
from WMAP3 + 2dF, dot-dashed from WMAP3 + SNLS and solid lines from the 
combination of all data. WMAP3-only constraints are very loose due to degeneracies 
between A and other parameters, especially fl m . 




Figure 2. Left: Joint 2D constraints in the Hq — A plane for a SUGRA potential. 
Lines references are the same as figure ^ Notice the strong degeneracy between the 
two parameters when using only CMB data. Marginalization over A yields a high 
value of f2 m with a large uncertainty. Adding additional data strongly suppress the 
high value of A and marginalized constraints are similar to results for ACDM. Middle: 
Joint constraints in the u> c — A plane. Right: Joint constraints in the as — A plane. 
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Table 2. RP parameters. 



Parameter 


WMAP3 


WMAP3 + 2dF 


WMAP3 + SNLS 


WMAP3 + 2dF 
+SNLS 


loofyi 2 

Q c h 2 
Ho 

T 

n s 

Log 10 (A DE ) 


9 9«O+0.083 

— U.Uoo 

n in9 c i+ - 0083 
59.71|? 

o.o9ii8;8ll 

0.96918:^ 
< 8.9 

u.ouy_ 103 

n 600+ 0080 

u.uuu_ 086 


9 OQO+0.072 

— U.U/4 

n nqqq+° o° 66 
u.uyyy_ o 6 5 

68.1+0 
0.09418:8} 5 

o.96ol8;81? 

< -3.6 
0.26618;8P 

0.67618:88^ 


9 ooc;+0.074 

— U.U (0 

1 01 ^+ - 0067 

U.lUiO„o.o067 

71.3+ 2 | 

0-09418:811 

0.958_ ; i 6 

< -7.9 
n 94c;+ -o 25 

n 7-1 O+0.057 
u - ' ±z, -0.053 


9 99^+0.070 
— U.U/ 1 

1 1 6+ - 0057 

U.lUlU_o 0056 

70.8+H 

o.o92i8:81 5 
0.95718:818 

< -7.7 
n 947+0.019 

u - z ^' -0.018 
n 7-1 O+0.050 
u - ' io -0.050 


Table 3. ACDM parameters. 


Parameter 


WMAP3 


WMAP3 + 2dF 


WMAP3 + SNLS 


WMAP3 + 2dF 
+SNLS 


lOOlVi 2 

n c h 2 

Ho 

T 

n s 

0-8 


9 oon+O- 073 
z,.z,z,u_ .072 

o in^4+ 0079 

U.1UO^_ 0080 

73.0+f 2 , 

o.o9oi8:8S 
0.95518:81^ 

949+0.034 
u - z ^ z -0.035 

7fi8+0-0 48 
u. / oo_ 048 


9 oon+O- 072 
z - zzu -0 .070 

o i n7o +aoo5 ° 

U.lU/U_o.0050 

72.2 + |? 

o.o86i8;81^ 

n osq+o.oi6 
u.yoo_ oi6 

0.24818:81 
o.767i8:8i? 


9 990+O.O72 
Z.ZZO_ .o72 

U.1U 1 O_ .0058 

72.2+11 

0.08718:811 

n q^+o-oie 
u.yoo_ oi6 

o.25oi8;8l 
o.77oi8:8 4 l 


9 91O+0.071 

z.ziy_ .070 

o.io77i8;8811 

71.91U 

o.o85i8;81^ 

n q^q+o-015 
u.yoo_ .oi6 
n 9C9+O.OI8 
u.z,dz,_ 018 
n 77n +0.036 
u - ' ' u -0.036 



this is not so for RP. We therefore confirm that SUGRA is still a viable cosmology; 
on the contrary, while RP is compatible with data, the cosmology it yields has just 
minor differences from ACDM. Let us then discuss first the SUGRA case, which leads 
to relevant shifts of some parameters in respect to ACDM. 

4.1. SUGRA potential 

A first point to be outlined is that there is no appreciable difference between the 
likelihood of SUGRA and ACDM, as the x 2 value, for SUGRA, is only very marginally 
improved. 

When only WMAP3 data are considered, the range of allowed Ade values is 
essentially unconstrained: A = log^A^/GeV) < 15, at the 95% confidence level. 
Considering the whole set of data, the limit is significantly more stringent: A < 2.1 
at 95% c.l. Assuming for the other parameters their expectation values, A = 2.1 
corresponds to a — 3.36. The strong difference from RP is due to the fact that, at 
late times, when approaches the Planck mass, the exponential term dominates the 
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Figure 3. Marginalized ID constraints on basic and derived parameters of RP 
cosmology, for the same combinations of data sets as in figure ^ 



15 
10 
5 

h 
-5 
- 10 



'1 1 1 1 1 1 1 1 1 1 1 1 








0.2 0.3 0.4 0.5 0.6 




0.08 0.1 0.1S 0.14 




Figure 4. Left: Joint 2D constraints in the f2 m — A plane for a RP potential. Lines 
references are the same as figure ^ 

potential, smoothing the dependence on the energy scale. 

As far as other parameters are concerned, mean values and confidence intervals are 
different from ACDM; shifts are more relevant for n s , Q m and Hq. The shift of the 
primeval spectral index n s , in particular, affects the impact of WMAP3 data on the 
nature of inflation (HI S3]- Using the whole data set, we find for ACDM: 

„ _ f| ncq +0.015 +0.031 /o\ 
ll, s — U.yjO_ 016-0 .030 i \°) 

while for SUGRA: 

n _ r, qfi o +0.017+0.038 (q \ 
IL S — U.aUO_g .018-0 .032 • \ u ) 

Here we have quoted the mean, the 68% and 95% c.l. errors; the 95% errors include 
the contribution from the 68% uncertainty. These figures show that a Zel'dovich 
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spectrum, excluded at the 95% c.l. (n s < 0.984), if ACDM models are considered, 
is compatible with the whole data set (n s < 1.001), at the same confidence level, for 
SUGRA cosmologies (see also |2"5]). 

More in detail, let us outline that, when using only WMAP3 data, constraints on 
Q m , H and as are very poor, in qualitative agreement with results for WMAP1 data [26J. 
In particular, H Q = 63 ± 7Km/s/Mpc, Q m = 0.32 ± 0.07 and a 8 = 0.63 ± 0.09. This is 
due to a strong degeneracy between A and fi m , which in turn affects estimates of H and 
as, as shown in figure El This degeneracy is qualitatively similar to the geometric w~Q m 
degeneracy in models with a constant equation of state. Fixing all other parameters, 
the effective value of (w) = J daw{a)VLD^,{a) / J daQD^ia), is an increasing function of 
Ade and the field becomes dominant at earlier times than in ACDM. When fitting CMB 
data, this behaviour yields higher fl m and lower as as Ade increases. Notice that there 
is no one-to-one relation Ade~w(ci), and results for models with constant w cannot be 
straightforwardly applied to constrain the energy scale. 

The main difference between WMAP1 and WMAP3 constraints is in the 
marginalized distribution of n s and r. The Integrated Sachs- Wolfe effect acts on low- 
i anisotropy multipoles more strongly in SUGRA than in ACDM. In WMAP1, this 
increases the r-n s degeneracy, and both r and n s are quite higher for SUGRA than for 
ACDM. WMAP3 estimates of TE and EE spectra allows a firm determination of r, and 
considering dynamical DE does not significantly worsen constraints on r and n s . 

Including other datasets in the analysis strongly reduces the allowed range for A DE , 
in turn significantly affecting the estimates of parameters degenerate with the energy 
scale. Both LSS and SNIa data exclude the highest values of A DE , corresponding to 
models with an high effective (w). The combination of WMAP3 and SNLS is particularly 
constraining for this class of models, and adding also 2dF data has only minor effects 
on estimates. In fact, current SN data seem to prefer models with equation of state 
around or slightly lower than -1, while for a quintessence model w(a) > — 1. 

SUGRA also favours higher values of u^, although differences are around 1/3 of 
the current uncertainty. As discussed above, Ade affects estimates of r and n s which 
in turn affect through the r — n s — u; b — A s degeneracy. However, if additional data 
allow to establish a clear difference in estimates for SUGRA and ACDM, Big Bang 
Nucleosynthesis measures could allow to distinguish the two models. 

Let us also outline that SUGRA models show an improvement in the goodness of fit, 
when LSS data are included in the analysis. Such improvement is modest, Ax 2 — — 1, 
and does not provide a compelling evidence in favour of SUGRA cosmology. However, 
the dynamical range of SUGRA models is still very wide and compatible with the energy 
scales of fundamental physics, and at present they offer a valid alternative to ACDM. 

4-2. RP potential 

Considering only WMAP3 data, results for RP are qualitatively similar to SUGRA. 
Upper limits on A DE are moderately tighter than in SUGRA, A < 8.9 (95% c.l.), and 
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discrepancies between RP and ACDM estimates are greater that those between SUGRA 
and ACDM. 

However, adding other data strongly suppress the high tail of the distribution on 
A DE ; using all data, we find A < —7.7 (95% c.L), well below the relevant range for 
EW and SUSY breaking. Moreover, estimates for other parameters are very similar to 
ACDM values, except for Q c h 2 and cxg. Both are lower than the corresponding ACDM 
figures by about 1 standard deviation. 

As already stated, previous works fHl EH! used various combinations of data to 
constrain the slope of RP potential, finding that a < 2. A direct comparison of the 
marginalized distribution on a and Ade is not immediate; however, we notice that 
assuming A DE = 10~ 7 GeV we find that 0.7 < a < 0.8 when f2 DE varies in the range 
0.1 — 0.9 and Hq = 70Km/s/Mpc. For these values of the parameters, the range of 
initial conditions for which reaches the tracking regime is very small, and RP models 
have fine tuning problems similar to ACDM |2Zj. Choosing initial conditions so that <j) 
does not reach the tracking regime by the present time, could alter the bounds on the 
potential parameters [2E]- However, without tracking, quintessence models lose one of 
their conceptual advantages over other DE candidates. 

Thus RP models, while still compatible with present data, do not seem to be a 
strong alternative to a standard ACDM cosmology. 

4-3. Summary 

We have used recent cosmological data to constrain the parameters of two relevant 
quintessence potentials, Ratra-Peebles (RP) and SUGRA. Both potentials depend 
on an energy scale, A DE , and the corresponding cosmologies reduce to ACDM for 
(A DE ) 4 ~ 10" 47 GeV 4 . 

Both for RP and SUGRA, current data allow to put only upper limits on A EE , due 
to a strong degeneracy between A DE and fl m . When considering only the third-year 
WMAP data, the degeneracy results in wider errors on several parameters, including Q m , 
H and <t 8 . When LSS and, especially, SN data are added to the analysis, constraints 
on Ad E strongly tighten. For RP potential, only values close to the cosmological limit 
are still allowed, and the model does not offer clear advantages to ACDM. 

For SUGRA potential, instead, values of A DE ~ 100 GeV, close to weak interaction 
energy scale, are still compatible with data. In addition, estimates of several parameters 
are different from the corresponding ACDM values. In particular, the spectral index, n s , 
is compatible with a Zel'dovich spectrum at less than 2 standard deviations. This could 
have a strong impact on constraints on the shape of the inflation potential. SUGRA 
models, thus, are still a viable and significant alternative to ACDM. 
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